NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



(NAS A-CR- 177 1 C9) TURBULENT BOUNDARY LAYERS N86-29149 

IK OSCILLATING FLOWS. PART 1: AN 
EXPERIMENTAL AND COMPUTATIONAL STUDY Final 
Technical Report, 1 Jun. 1982 - 31 Aug. 1985 Unclas 

(Iowa State Univ. of Science and Technology) G3/34 43241 


Iowa 


TURBULENT BOUNDARY LAYERS IN OSCILLATING FLOWS: 


AN EXPERIMENTAL AND COMPUTATIONAL STUDY 

William J. Cook 
February 10, 1986 


ISU-ERI-Ames-86434 

Final Technical Report (Part l'of 2 parts)* fc-; Grant NCC 2-200 

entitled "Study of Turbulent Boundary Layers in Oscillating Flows" 

Grant Period: June 1, 1982 to August 31, 1985 

William J. Cook, Principal Investigator 

J. D. Murphy, Grant Technical Monitor, NASA Ames Research 
Center, Moffett Field, CA 



*Part 2: Giddings, T. A. and Cook, W. J., "Response 
of Hot Element Flush Wall Gauges in 
Oscillating Flow," ISU-ERI-Ames-86435 , 
February, 1986. 


engineering 
research institute 

iowa State university College of Engineering, Ames, Iowa 50011 



i 


ABSTRACT 

An experimental-computational study of the behavior of turbulent 
boundary layers for oscillating air flovs over a plane surface with a small 
favorable mean pressure gradient is described. Experimental studies vere 
conducted for boundary layers generated on the test section vail of a 
facility that produces a flow vith a mean free stream velocity and a 
superposed nearly-pure sinusoidal component over a vide range of frequency. 
Flovs at a nominal mean free stream velocity of 50 m/s vere studied at 
atmospheric pressure and temperature at selected axial positions over a 2 m 
test length for frequencies ranging from 4 to 29 Hz. Quantitative 
experimental results are presented for unsteady velocity profiles and 
longitudinal turbulence levels obtained from hot vire anemometer 
measurements at three axial positions. Mean velocity profiles for 
oscillating flovs vere found to exhibit only small deviations from 
corresponding steady flov profiles, vhile amplitudes and phase 
relationships exhibited a strong dependence on axial position and 
frequency. Since sinusoidal flovs could be generated over a vide range of 
frequency, studies at fixed values of reduced frequency at different axial 
positions vere studied. Results shov that there is some utility in the use 
of reduced frequency to correlate unsteady velocity results. The 
turbulence level u' rms vas observed to vary essentially sinusoidally around 
values close to those* measured in steady flov. Hovever, the amplitude of 
oscillation and phase relations for turbulence level vere found to be 
strongly frequency dependent. Numerical predictions vere obtained using an 
unsteady boundary layer computational code and the Cebeci-Smith and Glushko 
turbulence models. Predicted quantities related to unsteady velocity 
profiles exhibit fair agreement vith experiment vhen the Cebeci-Smith 
turbulence model is used. 
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1. INTRODUCTION 

The behavior of turbulent boundary layers subjected to time-dependent 
free stream flows is of fundamental interest in several areas of fluid 
dynamics. Examples are those related to dynamic stall and buffeting that 
occur under certain conditions in the flow over helicopter rotor blades and 
blades in turboaachines . A number of experiments have been performed in 
recent years to study the behavior of turbulent boundary layers for flows 
with a periodic variation in free stream velocity. In parallel, 
theoreticians have worked to develop flow models and computational codes to 
predict the behavior of unsteady boundary layers. A number of these 
studies have been directed toward turbulent boundary layers developed on 
plane surfaces in essentially zero-mean-pressure-gradient flows as a first 
step in dealing with more complex flows. In the present investigation a 
joint experimental-computational study of the response of turbulent 
boundary layers for oscillating flows over a plane surface was performed. 
The purpose of this paper is first to describe the experimental study and 
the computational method and second to present comparisons of experimental 
results and predicted behavior for representative cases. Observations 
related to the comparisons are made. 

The concept of ensemble-averaging ('ielionis, 1981) provides the basis 
for experimental and theoretical treatment of the periodic turbulent 
boundary layer flows addressed here. The instantaneous value of a flow 
variable q at a fixed position in a periodic turbulent flow may be written 
as 

q(t) = q m + q(t) + q'(t) 

where q m is the time average or mean value of q, q(t) is the periodic 
component, and q'(t) is the turbulent fluctuation. The ensemble average of 
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q(t) is defined as 

q (0 - q m + q(t) 

vhere T 

- i l kw* 

q(t) is a function defined by J values of q each of which is given by 

i N 

q< t ± ) - £ 5 ql t A + (n-i)yj 
n-1 

In the above equation, y is the period of the free stream oscillation, N is 
the number of cycles over which the average is obtained, and tj - (i/J)y, i 

* 1 f ■ > * , J* 

The ensemble-averaged description of the free stream velocity in an 
oscillating flow with a sinusoidal variation for U(t) is, at a fixed 
position in the flow direction, 

U(t) - U m + UjCos «t (1.1) 

and ensemble-averaged velocity in the boundary layer is 

u(y,t) - u m (y) + G(y,t) = u m (y) + u 1 (y)CosIwt + $(y)] (1.2) 

where the subscript 1 denotes the half amplitude of the velocity 
oscillation, o> is the angular frequency, and ♦ is a phase angle. The ratio 
u^/U^ may exceed unity in the boundary layer and $ may exhibit both 
positive and negative values. 

2. EXPERIMENTAL FACILITY 

2.1 Facility Description 

A number of experimental facilities have been used to study unsteady 
turbulent boundary layers in oscillating flows with zero or essentially 
;:ero mean pressure gradients. A review of facilities used in oscillating 
flow studies has been prepared by Carr (1981). Some of the facilities used 
by other investigators with air as the test medium in zero-mean-pressure 
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gradient flows will be briefly discussed. The first significant study of 
turbulent boundary layers is oscillating flows was performed by Karlsson 
(1959), and was carried out using an oscillating flow facility consisting 
of an upstream blower, screens, a contraction section, and a test section 
of rectangular cross section. The boundary layer generated on the test 
section wall was studied. Rotating vanes installed near the test section 
exit generated oscillating flow in the test section. Variations of this 
flow-generating technique have been used by others including Simpson, et al 
(1981), and Cousteix, et al (1977). Other successful methods for 
generating oscillating flows in air have been employed by Brembati (1975), 
and Patel (1977). Brembati used a tunnel with a flexible test section wall 
that oscillated sinusoidally normal to the flow direction to create a 
variation in free stream velocity. Patel used a blower-driven contraction 
section that delivered flow to a horizontal test plate enclosed only by 
vertical walls. Flaps at the exit of the flow-delivery section were 
oscillated sinusoidally to induce travelling vortices in the flow over the 
test plate, thereby producing a unique type of oscillating flow. With the 
exception of Patel's facility, all of the above mentioned facilities 
produce oscillating flows that are relatively free of higher harmonics only 
when operated at a discrete frequency or, in some cases, at approximately 
twice the discrete frequency. 

The facility developed in this study generates flows by a method that 
is conceptually quite different from those used in other studies. The 
facility is illustrated in Figure 1 and is further described in Table 1. 

It consists of an entrance section and a test section of square cross 
section and length L which is terminated by a convergent-divergent nozzle 
that has a two-dimensional contour. The nozzle discharge region pressure 
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is sufficiently low during testing to produce sonic flow at the nozzle 
throat. When the position of the wedge of the wave generator is fixed , the 
test section flow, which is subsonic, will be steady. Vhen the wave 
generator mechanism is operated, the test section free stream flow at any x 
position is characterized by a non-zero mean velocity and a superposed 
oscillating component. The mechanism consists of a scotch yoke which 
imparts sinusoidal axial motion to the wedge when the input shaft is driven 
at constant speed. As the wedge undergoes motion, the sonic throat area 
varies, thereby varying the mass rate of flow and producing a wave pattern 
in the test section. The wave form developed depends on the contour of the 
wedge surface. A ramp-type symmetric wedge with a small leading-edge 
included angle produces an essentially pure sinusoidal velocity variation. 
For a given nozzle-wedge configuration, a change in wedge stroke produces a 
change in the amplitude of free stream velocity. Vedge strokes used during 
the course of this study are listed in Table 1. 

The turbulent boundary layer studied in the present investigation was 
that generated on the test section wall. Boundary layer trip wires placed 
on all four walls at x * 0.13 m (Figure 1) were chosen to cause transition 
to occur at the trip location. The wire diameter to produce transition at 
the trip was determined by the method of Smith and Clutter (1959). The 
test section was square in cross section (7.5 cm x 7.5 cm) for its full 
length. As a result of boundary layer growth, a mild favorable pressure 
gradient existed in the flow. The origin of the wall boundary layer was 
taken at the trip location. The nominal mean test section velocity was 50 
m/s. Thus the flow was in the incompressible regime. 

The present technique for generating oscillating flows has the 
following advantages. First, the sonic throat isolates the test section 
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from dovnstream disturbances. A properly designed entrance section will 
deliver flow to the test section with a relatively lov disturbance level. 
The result is a test section flow that has a minimum of extraneous 
disturbances. Second* with a ramp type wedge an essentially-pure 
sinusoidal free stream velocity variation can be produced over a vide range 
frequency. Third, independent changes in the mean flow velocity and the 
amplitude of oscillation can be made with relative ease by changing the 
configuration of the nozzle throat and the wedge of the wave generator. 
Flows with large amplitudes of oscillation can be generated. Fourth, due to 
the manner in which the flow is generated, an analytical description of the 
mainstream flow can be obtained. An accurate description of the mainstream 
velocity distribution with position and time is required for theoretical 
analysis of the boundary layer behavior. The following is the development 
of this description. 

2.2 Analytical Description of Free Stream Flow 
The oscillating free stream flow is assumed to be isentropic and 
effects of the wall boundary layer are ignored. Thus the flow is described 
in terms of the space coordinate x (Figure 1) and time. Under these 
conditions the free stream velocity in the test section is 

U(x,t) - U m + U(x,t) 

The analysis presented here is for small amplitudes of oscillation. 

Therefore, the amplitude of the varying component U(x,t) is assumed to be 

small compared tu the mean value U m> The governing equation for the 

velocity variation that develops in the test section under these conditions 

is the wave equation in the form 

(a 2 - U 2 )$ - 2U * „ - = 0 (2.1) 

' m m xx m xt tt ' 

where t is the velocity potential function, which is related to U(x,t) by 
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0(x,t) ■ 3f(x,t)/3x 

Solutions to equation (2.1) have been obtained in this study for sinusoidal 
U(x,t). For reasons related to boundary conditions for equation (2.1) the 
solution vill be written in terms of the Mach number M(x,t) - U(x,t)/a m 

where a is the mean acoustic velocity. The Mach number in the test 

lu 

section is 


M(x,t) - M m + M(x, t) 


» M m - S(x)Cos[wt + A(x)J (2.2) 

where to is the wedge angular frequency, S(x) is the half amplitude of the 
Mach number variation, and A(x) is a phase angle. The latter are given by 

(2.2a) 

(2.2b) 


9 2 1/2 

S(x) - [ (BCosK + CCosnKr + (BSinK - CSinnKrl '* 


»/ v \ a -1, BSinK - CSinnK , 

A(x) - tan l BCosK + C cosnK^ 


where 

n - (1 - M m )/(1 + M m ) (2.2c) 

K * (xw/a m )/(l - M m ) (2. 2d) 

and B and C are constants yet to be determined. The corresponding pressure 
in the test section is 


P(x,t) . P m + P(x,t) 

where P(x,t) * p(x,t)/kp m is a dimensionless pressure. p(x,t) * P m + 
p(x,t) is the absolute pressure, k is the specific heat ratio yc v tor 
the test gas (air in this study) and p m is the mean pressure. The solution 
for the pressure is 


in which 


and 


P(x,t) = P m + W(x)Cos[wt + <Kx)] (2.3) 

2 2 1/2 

W(x) = [ (BCosK - CCosnK) +(BsinK + CSinnK) ]' (2.3a) 
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««> - tan ' 1 Sco^k - SS3 < 2 - 3b > 

For a specified value of the mean Mach number, equation <2. 2c) yields n in 
equations (2.2a,b) and (2.3a,b). Additionally, equation (2. 2d) relates K 
to x when wand a m are specified. Equations <2. 2) and <2. 3) then yield the 
time and position dependent Mach number (velocity) and dimensionless 
pressure for the test section provided the constants B and C in equation 
(2.2a,b) and (2.3a,b) can be determined. Evaluation of these constants 
requires two boundary conditions. Assumptions related to these boundary 
conditions are that the entrance section flow and the flow in the 
converging section of the nozzle are isentropic and quasi-steady. These 
assumptions are appropriate if the entrance section and the converging 
section of the nozzle are short compared to the length of the test section. 
The first boundary condition is obtained at x ■ L. The ratio of the cross 
sectional area at the nozzle entrance, A^, to the sonic are A * changes with 

•ff 

time as A changes due to the motion of the wedge. Imposing the assumption 
of isentropic quasi-steady flow in the converging section of the nozzle, 
M(L,t) is given implicitly by the relation between A/A and the Mach number 
for isentropic flow, which is 

A L /A*(t) - [l/M(L,t)]{[2/(k-l)] + 

[(k-l)/(k+l)HM(L,t)] 2 )^ k+1)/l2<k " 1)1 (2.4) 

For small changes in A at low subsonic values of Mach number, M is nearly 

jlf ^ 

linear with A /A. Therefore, a sinusoidal variation in A for low Mach 
number flows will produce a variation in M(L,t) that is essentially 
sinusoidal. From equation (2.2) evaluated at x = L, 

M(L, t) = M m -S(L)Cos[wt + A(L) ] 

where S(L) = represents the half amplitude of the Mach number variation 
at x = L. M m is found from the A/A - Mach number relation using A/A 
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determined with the wedge positioned at mid-stroke. Further, depends 
only on the geometry of the wave generator, i.e,, ths nozzle throat cross 
sectional area, the wedge included angle, and the wedge stroke. Thus, for 
a specified configuration of the wave generator, the boundary condition at 
x - L yields one equation relating the constants B and C. This equation is 
obtained from equation <2. 2a) as 

S L - KBCosKl + CCosnl^) 2 + (Bsin« L - CSinni^) 2 l 1/2 (2.5) 
The second boundary condition is related to the pressure variation at 
x * 0. Equations (2.3a) and (2.3b) evaluated at x ■ 0 yield 

W(0) - V Q - B - C (2.6) 

'KO) « 0 

From equation (2.3) the pressure variation at x ■ 0 becomes 

P(<M> « W Q Cos cot (2.7) 


Equation (2.6) provides the second relation between the constants B 
and C. From equations (2.5) and (2.6) the following expressions are 
obtained 

B - [W 0 + (W Q 2 - 2<V o 2 - S L 2 )/(l + z)} 1/2 ]/2 

C - B - W Q (2.8) 

where 


z ■ CosKjXosnKk - Sin^SinnK^ 

While S L is known from specification of the wave generator geometry, W Q is 
unknown. Obtaining the value for W Q is related to the quasi-steady 
isentropic flow assumption for the entrance section. This assumption 
relates p(0,t) to M(0,t) through the relation 

p(0,t)/p° = {l/[l+(k-l)(M 2 (0,t)/2)]) k/(k " 1) (2.9) 

where p° is the stagnation pressure of the air in the supply reservoir; the 
atmosphere in the present case. The expression 
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p(0,t)/p° - p m /p° + p(0,t)/p° 

combined with equation (3.9) tt. form the dimensionless pressure variation 
at x « 0 yields 


*<0,t) - p(0, t)/kp 


m 


(p°/kp„) 


[(: 


l+(k-l) (M 2 (0, t)/2' 


k/(k-l) 




( 2 . 10 ) 


Equation (2.10) yields a periodic but not exactly sinusoidal P(0,t) for a 
sinusoidal H(0»t). Denoting P' as the pressure variation obtained from 
equation (2.10), an equation similar in form to equation (2.7) can be 
written as 

P'(0,t) - V 0 'Cos«t (2.11) 

where W ' is the half amplitude of the pressure variation produced by 
equation (2.10). M(0,t) for equation (2.10) is given by equation (2.2) as 

M(0, t) »• M m + M(0,t) » M m - S(0)Cos[u>t +A(0) ] 

- M - S Coswt (2.12) 

mo 

However , evaluation of S Q from equation (2.2a) requires values for B and C, 
which according to equations (2.8) depend on the known and the unknown 
tf Q . The solution is obtained through iteration by first guessing a value for 
V Q . This will yield values for B and C through equations (2.8). S Q is then 
found from equation (2.2a); S Q - B + C. Then substitution of M(0 f t) from 
equation (2.12) into equation (2.10) will then yield V Q ' in equation (2.11), 
which is then* compared with the guessed value of U Q . Iteration by means of 
a computer will result in V Q ' * W Q within a specified small tolerance, at 
which point values for B and C are established. These values then produce 
the analytical predictions for the Mach number (velocity) and pressure by 
means of equations (2.2) and (2.3). Comparisons of predicted and 
experimentally observed free stream flows are made subsequently. 
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3. DATA ACQUISITION SYSTEH 

A hot vire anemometer system was used as the primary means of 
obtaining data. Figure 2 describes the data acquisition system used 
throughout this study. A DISA P14 single-wire probe mounted on a 4 mm 
diameter probe support was positioned through a vail of the test section. 
All data vera taken at the half-height position on a vertical vail. The 
hot vire system vas calibrated using a TSI model 1125 airflow calibrator, 
and the calibration vas checked periodically during the experiments. 

During the course of this investigation several different hot vire probes 
were calibrated and used. The data acquisition system vas used for both 
steady and unsteady flow measurements. For steady flow boundary surveys 
the system vas programmed to read and process 500 hot wire voltage values 
at each probe position. For oscillating flow studies the system vas 
programmed to read 20 hot vire voltage values at equally-spaced time 
intervals per cycle of flov oscillation and vas triggered electronically 
once per revolution of the wave generator drive shaft. The shaft frequency 
vas monitored by a photoelectric device. The system vas programmed to take 
data for 210 cycles. Ensemble averaged velocities were obtained by 
processing the 210 hot vire voltage readings at each of the 20 points in 
the cycle. The data processing included the statistical fitting of a 
cosine wave to the 20 ensemble averaged values. The maximum data 
acquisition rate for the system vas 1000 voltage readings per second. Data 
taken during this study were acquired at rates veil belov this value. 

4. STEADY FLOW EXPERIMENTS 

Detailed hot wire measurements were made at various x stations along 
the test section vith the wave generator wedge locked at mid-stroke to 
provide a description of the test section flow under steady conditions. 
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Figure 3 presents as a function of test section axial position results for 
several flov quantities determined from the steady flow measurements. Due 
to the growth of the wall boundary layers there was a detectable increase 
with x of the steady centerline velocity U . The velocity increase 
described in Figure 3 corresponds to a measured pressure drop of 4.1 cm of 
water (0.40 kPa) over the test section length. Thus a mild favorable 
pressure gradient existed in the test section. Figure 3 also shows the 
test section centerline turbulence intensity u' r /U s vs x, where u' r - u' rmg 
. [u'2]'* . These measurements were made by using a true rms meter and a 10 
kHz low-pass filter in place of the voltmeter and the computer in the data 
acquisition system in Figure 2. The purpose of the filter was to eliminate 
high frequency electrical noise on the anemometer output signal that 
produced false readings at low turbulence levels. The turbulence intensity 
ranged from 0.1 percent at the test section entrance to 0.4 percent at x « 
1.7 m. 

Measured values for 8 , the boundary layer velocity thickness for 

s 

steady flow determined at 99 percent of the local centerline velocity, are 
also shown in Figure 3. The curve fitted through the data points indicates 
that the boundary layer grew according to the relation S g - 1.475 ’ » 

where 8 is in cm and 5 is the distance in meters from the boundary layer 

s 

trip wires (Figure 1). The boundary layer growth expression is consistent 
with the presence of a mild favorable pressure gradient. For zero pressure 
gradient turbulent flows, for which the accepted exponent is 0.8, a more 
rapid growth with 5 occurs. 

The ratio of 8 to the test section height h increases with 5 to a 
s 

value of about 1/4 at the most down-stream measurement station, 5 = 1-58 m. 
The thickness of the wall boundary layer relative to the test section 
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height h is of concern in rectangular and square ducts since for 
sufficiently large values of S/h secondary flows in the corners apparently 
cause cross flow in the boundary layer which destroys the nominally two 
dimensional boundary layer flow at mid plane. This matter has been studied 
experimentally and analytically by Brederode and Bradshaw (1978). It was 
found that the effect of secondary corner flows on the mid-plane boundary 
layer is negligible for 8/h < 1/4, and that under this condition, the 
boundary layer flow on the wall at mid-plane in a square duct of constant 
cross section is quite close to that on a flat plate. Since the condition 
8/h < 1/4 is nominally met in the present study, flat-plate-like flows were 
expected on the test section side wall for £ < 1.58 m. 

Momentum thickness Reynolds numbers shown in Figure 3 for the boundary 
layer exceed 3600 for x > 0.74 m and indicate that a fully turbulent 
boundary layer was generated. The values shown for the boundary layer 
shape parameter H = 8*/0 are typical of those for a flat plate turbulent 
boundary layer in a zero pressure gradient flow. A small pressure gradient 
does not significantly influence H. 

Figure 4 presents measured steady flow velocity profiles and 

turbulence intensity in terms of y/8 for three 5 stations; 0.61, 1.12, and 

s 

1.58 m. 8 was taken as that given by the equation fitted to the 
s 

experimental values of 8 vs £ in Figure 3. Figure 4a shows velocity 

o 

profiles in terms of u/U and y/8 . The three profiles shown are in good 

s s 

agreement. Figure 4b shows steady boundary layer velocity profiles at the 
same three £ stations in terms u + vs y + where u + = u/u T , u T = (%/p) > 

and y + = u^y/v. t w was evaluated by the Ludwieg-Tillman law (Coles, 1968). 
The experimental results at 1.12 m and 1.58 m agreement at low values of y + 
with the logarithmic law region equation. A small difference exists 
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bet ween the experimental and predicted results at £ - 0.61 m. 

Boundary layer turbulence intensities in terms of u' /u and y/6 are 
presented in Figure 4c for the three £ stations. Turbulence intensities in 
the figure were determined by computer processing of the 500 hot-wire 
voltage values read using the data acquisition system in Figure 2 and, 
except at low values of turbulence intensity, were in agreement with 
readings obtained using the previously-mentioned 10 kHz low pass filter and 
rms voltmeter. As expected, the results at the three values of £ collapse 
toward a single curve. 

The collective results in Figures 3 and 4 indicate that a turbulent 
boundary layer typical of that developed in a flow with a mild favorable 
pressure gradient was produced in the test section under steady flow 
conditions. 

5. OSCILLATING FLOW EXPERIMENTS 

5.1 Free Stream Velocity Variation 

As noted previously, the facility used in the present study produces 
an essentially-pure sinusoidal velocity variation in test section for a 
wide range of frequency. This is demonstrated in Figure 5 which shows 
velocity variations measured in the test section for two values of 
frequency at two positions using the data acquisition system described in 
Figure 2. The wave generator configuration described in Table 1 with a 
wedge stroke of 1.91 cm was used to produce the velocity variations shown 
in Figure 5. The symbols represent velocity values obtained by ensemble 
averaging velocity measurements (20 points per cycle) over 210 cycles. The 
curves are cosine waves fitted to the 20 ensemble averaged velocity values 
for the cycle. The three upper curves describe velocities measured in the 
free stream at two positions, £ = 0.305 and 1.58 m, while the lower curve 
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1s the velocity measured in the boundary layer at y/5 - 0.1 at 1.58 m. It 

s 

is evident in each case that the curve closely £its the data. Expressions 
for the fitted curves in the form of equations (1.1) and (1.2) permitted 
accurate evaluation of the mean velocities, amplitudes, and phase angles 
from the experimental data. 

The two upper curves in Figure 5 demonstrate the effect of frequency 
on the free stream velocity variation at a fixed value of £. These curves 
are for velocity variations at 5 * 1.58 m. It is evident that the 
amplitude for the velocity at 29.4 Hz is larger than that for 15.5 Hz. 
Further, the velocity variation at 29.4 Hz leads that at 15.5 Hz by a small 
angle. The tvo free stream velocity curves in Figure 5 obtained at 15.5 
Hz, one at 5 * 0.305 m and the other at 5 - 1.58 m, demonstrate the effect 
of position on the free stream velocity at a fixed frequency. Due to the 
increase in the mean free stream velocity in the direction of flow (as 
noted also for steady flow in Figure 3), the curve for £ « 1.58 m at 15.5 
Hz lies above that for £ = 0.305 m in Figure 5. However, at £ * 1-58 m 
is smaller than at 0.305 m. It is also evident from the figure that the 
free stream velocity at 1.58 m leads the velocity at 0.305 by a small 
angle. 

The variation in the free stream velocity amplitude and phase angle 

with frequency and position can be predicted for a specified configuration 

of the wave generator by the analytical method previously outlined. For 

this analysis it is necessary to specify the mean Mach number M m and the 

Mach number amplitude at x = L, S. in equation (2.5). M can be found from 

l m 

equation (2.4) using A(L)/A determined from the wave generator 
configuration with the wedge positioned at mid-stroke. The extremes of 
A(L)/A* can be determined from the wave generator configuration for a 
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speci£ied wedge stroke, from which can then be found using equation 
(2.4). Solutions for S(x) and A(x), equations (2.2s) and (2.2b), can then 
obtained at the mean Mach number M m . It may be difficult to accurately 
determine M m and from the wave generator configuration if small 
dimensions are involved. An alternate method for evaluating M m and 
involves measuring the free stream velocity in the test section with the 
wedge fixed at each of three positions; fully upstream, mid stroke, and 
fully downstream. From these measurements M m and S L can be accurately 
determined. Figure 6 shows solutions for S(x) and A(x) at four values of 
frequency ranging from 1 to 29.4 Hz for the wave generator described in 
Table 1 at a wedge stroke of 1.91 cm. The values of M m and for the 
results in Figure 6 are 0.142 and 0.018 and were obtained using the 
alternate method from velocity measurements made at the test section 
half-length position. Also shown in Figure 6 are values of S and A 
obtained from velocity measurements in the free stream at three values of 
frequency at several x positions in the test section, i.e., from results 
like those for the free stream velocity in Figure 5. The agreement between 
the analytical solution and the experimental results is generally good. 

Both the experimental and analytical results indicate that the amplitude of 
oscillation departs from low frequency uniformity with x as frequency 
increases and that the phase angle exhibits a progressively increasing 
dependence on x as frequency increases. The results in Figure 6 
demonstrate that both the free stream velocity variation and the phase 
angle variation with frequency and x can be predicted with reasonable 
accuracy by the analytical method. 

5.2 Boundary Layer Velocity Variation 
The major emphasis of the experimental phase of this investigation was 
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on acquiring data to describe the boundary layer behavior in oscillating 
flows. A total of twenty boundary layer surveys were conducted at four 
test section axial positions, 5 - 0.61, 0.91, 1.12, and 1.58 m, for 
frequencies ranging from 3.85 to 29.4 Hz using the data acquisition system 
described in Figure 2. Frequencies in this range are considerably below 
the estimated turbulent burst frequency, 500 Hz, obtained from the 
expression U m /58 s (Rao, Naradinha, and Badri Narayanan, 1971). The results 
of each survey are denoted as a data set. Data set designations and 
conditions at which each was obtained are tabulated in Table 2. Conditions 
for the surveys were selected to provide results for comparisons at fixed 
values of position, frequency, and reduced frequency 5> ■ w£/U m< The range 
of reduced frequency covered by the surveys was 0.54 to 4.16. Vhile u is a 
similarity variable for zero pressure gradient laminar oscillating boundary 
layers, exact similarity in terms of w for turbulent oscillating boundary 
layers does not exist (Telionis, 1981). Nonetheless, data were obtained at 
fixed values of <a at several $ stations in an attempt to assess the extent 
to which oi serves as a correlating variable. 

The amplitude of oscillation of the free stream velocity may have an 

influence on boundary layer behavior. It is evident from Figure 6 that the 

amplitude of velocity oscillation at high frequencies varies with position 

along the test section. The values for U./U listed in Table 2 were 

1 m 

obtained from free stream velocity measurements at the location where the 
survey was conducted. U^/U m could be varied from survey to survey by 
changing the wedge stroke. The wedge stroke for each survey is listed in 
Table 2. There are four pairs of data sets for which conditions were 
identical or essentially identical within the pairs except for wedge 
stroke. These pairs of surveys were performed to assess the effect of 
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amplitude of velocity oscillation on results for a limited range of Uj/U m< 
Results for the oscillating flow boundary layer surveys vill be 
presented in terms of the quantities in equations (1.1) and (1.2). 

Equation (1.2) represents the ensemble averaged velocity in the boundary 
layer produced by the free stream velocity variation described in equation 
(1.1). Descriptive quantities are u m /U ra , the ratio of the mean boundary 
layer velocity to that in the free stream; u^/Up the ratio of the 
half-amplitude of oscillation in the boundary layer to that in the free 
stream; and the phase angle of the boundary layer velocity variation 
relative to that in the free stream. It is useful to present these 
quantities in terms of y/8 , where 6 is the local boundary layer velocity 
thickness for steady flow as presented in Figure 3. Use of this ratio 
readily permits evaluation of the extent of unsteady effects in the y 
direction relative to the steady boundary layer thickness. Values for U m 
and in the free stream and u m , Up and ♦ in the boundary layer were 
determined from equations for the cosine curves fitted to ensemble averaged 
velocity data. A typical case is presented in Figure 5, in which the 
boundary layer velocity variation at y/S » 0.10 and the corresponding free 

5 

stream velocity variation at £ * 1.58 m and 15.5 Hz are shown. It is seen 
from the figure that the velocity variation in the boundary layer is 
essentially sinusoidal and that it leads the velocity variation in the free 
stream. 

' Figure 7 presents representative results for mean velocity profiles in 

terms of u /U, and y/S at three values of £. The oscillating flow data 
mm's 

cover a wide range of frequency (see Table 2). Results of the steady flow 
velocity surveys are also shown in Figure 7. Estimates of uncertainties 
related to the velocity measurements indicate a small uncertainty region 
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about the size of the symbols. At ( » 0.61 m, Figure 7a, the mean profiles 
for oscillating flow agree yell with the steady flow profile. Results at 5 
* 1.12 and 1.58 m, Figures 7b and 7c, show small differences in the 
unsteady and steady mean profiles, with the difference in Figure 7c the 
most noticeable. 

Figures 8 and 9 present results for the amplitude ratio u^/U^ and the 
phase angle f respectively at £ - 0.61, 1.12, and 1.58 m in terms of y/6 

s 

for several values of frequency of oscillation. The results included in 
these figures were obtained for flows generated with the wedge stroke set 
at 1.91 cm, and, as noted in Table 2, values of U^/U m ranged from 0.121 to 
0.239. Prior to discussion of these figures it is appropriate to consider 
the uncertainties of the results presented therein since these bear on 
conclusions drawn. 

One of the sources of uncertainty in the oscillating flow results is 
related to the number of cycles over which the velocity measurements are 
ensemble averaged. Generally, due to the random nature of turbulence, 
ensemble averaging over a small number of cycles results in larger 
uncertainties for u^ and ♦ than does averaging over a large number of 
cycles. In order to assess these uncertainties, the input to the data 
acquisition system (instantaneous readings of the oscillating flow 
velocity) was modeled using a digital computer and a statistical analysis 
was performed at various positions in the boundary layer for oscillating 
flows with ensemble averaged velocity variations, as described by equation 
(1.2), typical of those observed experimentally. The level of turoulence 
is an important parameter in the uncertainty analysis. In the free stream 
where u' r , the rms turbulence level, is small compared to the amplitude of 
velocity oscillation U^, an accurate representation of the ensemble 
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averaged velocity variation can be obtained by averaging over a relatively 
few cycles, while in the boundary layer, where u' r may be near u^, data £or 
a much larger number of cycles must be acquired for accurate evaluation of 
the ensemble averaged velocity variation. In the uncertainty analysis the 
flow was modeled for 210 cycles, the number of cycles over which 
experimental data were ensemble averaged. A Gaussian distribution about u 
was assumed for u' and u' r at a given y/8 s was taken as that measured for 
oscillating flows. (The variation in turbulence level measured for 
oscillating flows is discussed in a later section.) Next, at each of the 
20 points in the cycle, 210 u' values were generated by random entries onto 
the Gaussian distribution curve defined by the value of u' . These were 
superposed on the assumed ensemble averaged velocity variation, thus 
producing data that simulate the instantaneous velocity values determined 
from hot wire measurements. These data, when ensemble averaged, yield 
values of u m , Uj, and t, which differ from the input values because of the 
superposed random u' values. By repeating this simulation process at a 

given y/5 a large number of times, the uncertainties in terms of standard 

s 

deviations for u m , u^ and + could be estimated. Figure 10 presents 
uncertainties determined in this manner for u^ and $ in terms of y/& s for 
flows produced with an oscillator wedge stroke of 1.91 cm. The uncertainty 

for u. is not shown since it is an order of magnitude smaller than that for 

m 

u^. As expected, in the free stream where the turbulence level is small 
compared to u^, the uncertainty is small for both u^ and 4 >. As the wall is 
approached, u' r increases and u^ decreases, resulting in larger 
uncertainties in u^ and $. 

Other sources of uncertainty for the results for oscillating flows 
include uncertainties in frequency of oscillation, & , hot wire position y, 

O 
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and, for u^/Up the hot wire calibration. Use of a dial micrometer for 

positioning the hot wire resulted in very small uncertainties in y. The 

uncertainty in 5 vas estimated to be ±4 percent. Thus, the uncertainty in 

y/S vas ±4 percent. Estimates of the uncertainty of results were made by 

combining the uncertainties from the sources using the method of Kline and 

McClintocIc (1953). The uncertainty range (+a) in Figure 8 for u^/U^ from 

all factors vas estimated to be the symbol size (±2 percent) at large y/5 

and tvice the symbol size (±4 percent) at small values of y/8 . Hence, the 

s 

uncertainty for Uj/U^ is a relatively small. 

The results for in Figure 9 involve larger uncertainties. An 
analysis of the uncertainty in 4 due to an uncertainty of ±0.5 percent in 
frequency of oscillation for the experiments yielded ■ ±0.5 degrees. 
Combining of this uncertainty vith that from Figure 10 for $ yielded an 
estimate of the total a ^ that ranged from ±0.5 degrees in the free stream 
to ±1 degree at y/5 * 0.05. Bars indicating the uncertainty range on $ 
are shovn on selected curves dravn through the data points in Figure 9. 

Results in the form presented in Figures 8 and 9 permit assessment of 
the boundary layer behavior at fixed axial positions as frequency varies. 
Several observations can be made from these results. The results for u^/U^ 
vill be considered first. It is noted in Figure 8 that the influence of 
unsteadiness extends beyond the steady boundary layer thickness to values 
of about 1.48 g at each 5 position. Significant changes in u^/U^ profiles 
occur at a fixed position as frequency changes, as evidenced for example by 
the profiles at £ = 1.12 m, Figure 8b. As frequency increases, u^/U^ is 
larger than unity in an increasing portion of the boundary layer and the 
maximum value of u^/U^ increases. This trend continues to a frequency of 
21.3 Hz. The profile at the largest value of frequency, 29.4 Hz, exhibits 
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a smaller maximum value of u^/U^ that occurs nearer the vail. The 

collective results in Figure 8 can also be discussed in terms of the 

reduced frequency. At low reduced frequency (data set A1 in Figure 8a and 

data set C2 in Figure 8b, for which 5> ■ 0.54) the profiles exhibit u^/U^ 

values that are slightly larger than unity in a small portion of the 

boundary layer. Beyond Q ■ 0.54 for flows with <3 ranging to 3, the 

profiles show u^/U^ > 1 for increasingly larger portions of the boundary 

layer. There is also a trend to larger maximum values for Uj/U^ that occur 

at smaller values of y/&_, For flows with w> 3 (data set C6 in Figure 8b 

s 

and data set D6, Figure 8c) the fraction of the boundary layer for which 
u^U, > 1 decreases, as does the maximum value of Uj/U^ and the value of 
y/& s at which the maximum occurs. 

The profiles for the phase angle 4, Figure 9, show an influence of 
frequency of oscillation to y values of about 1.56 and exhibit pronounced 

S 

effects of frequency. The results for £ ■ 1.12 m in Figure 9b are somewhat 
typical. At the two lowest frequencies 4 is positive throughout the 
boundary layer and exhibits a maximum at about y/6 - 0.3. As frequency 

S 

increases, regions of negative 4 develop toward the outer part of the 
boundary layer. Quite large negative values occur at the higher 
frequencies near y/8 * 0.8 As y/8 decreases, 4 becomes positive, and at 
a frequency of 21.3 Hz, the curve clearly exhibits a maximum value at a low 
value of y/6 . Due to the uncertainty in 4 of ±1 degree at low y/8 and 
the sparseness of the data near the wall, the tendency to exhibits maximum 
is less clear for the results at a frequency of 29.4 Hz (data set C6). The 
results at £ * 1.58 m, Figure 9c, exhibit characteristics similar to those 
at ^ * 1.12 m, while at $ = 0.61 m (Figure 9a), the profiles generally show 
only increasing values of 4 as the wall is approached. 
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The results in Figures 8 and 9 were obtained for flows generated with 

a wedge stroke of 1.91 cm. Table 2 lists four data sets that were obtained 

with shorter wedge strokes, and therefore had smaller amplitudes of 

oscillation. Except for amplitude of oscillation, the parameters for these 

four data sets, AS, C5, D5, and D7, respectively match or closely match 

data sets A4, C4, D4, and D6» Comparisons of the matched data sets are 

presented in Figures 11 and 12. The profiles of UjVUj for data sets A4 and 

A5 for which £ » 0.61 m and the frequency is 21 Hz are shown in Figure 11a. 

Although the local value of Uj/U m for data set A5 is two-thirds of that for 

data set A4, there is only a small difference in the profiles and this 

occurs at low values of y/6 . U.,/U m for data set C5 is about two-thirds of 

that for data set C4. The comparison in Figure lib for data sets C4 and C5 

show no detectable influence of amplitude of oscillation. The profiles for 

data sets D4 and D5, Figure 11c, are also essentially identical. The local 

value of Uj/U m for data set D5 is 60 percent of Uj/U m for data set D4. 

Figure 11c also shows profiles for data sets 06 and D7. These profiles, 

for which the frequency is 21 Hz, exhibit a difference for the full range 

of y/6 . The local value of U./U for data set D7 is one-half of that for 
s l m 

data set D6. 

Comparisons showing the influence of amplitude of oscillation on $ 
profiles are presented in Figure 12. The profiles for data sets A4 and A5 
show a small but distinct difference near y/5 g ** 1. Near the wall, the 
difference is about the size of the uncertainty in the angle measurements. 
(Due to the smaller amplitude of oscillation relative the the turbulence 
level in the boundary layer, uncertainties in results obtained at the 
shorter wedge strokes are larger than those for the 1.91 cm stroke.) The 
same situation exists in the comparisons of data sets C4 and C5, Figure 
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12b, and D4 and D5, Figure 12c. Within the uncertainty o£ the results, 
the profiles for data sets D6 and D7 show little difference over the full 
range of y/&_. Inspection of Figures 12 and 13 leads to the conclusion 
that for the flows investigated, the amplitude of velocity oscillation has 
at most a small influence on the boundary layer behavior. 

Figures 13 and 14 present profiles of Uj/U^ and f in terms of y/S $ for 
three values of frequency; 7 Hz, 15 Hz, and 21 Hz. All data sets at or 
near these frequencies are included in the figures. For data sets where 
conditions were the same except for wedge stroke and for which results were 
nearly the same, a single curve has been drawn through the collective data 
points for clarity. The results in Figures 13 and 14 allow visualization 
of the evolution of profiles for u^ and 4 ms functions of axial location at 
fixed values of frequency. In interpreting these results it is necessary 
to take into account the phase shift in ■* free stream velocity that 
occurs in the direction of flow at a fixed frequency as described in Figure 
6 . 

As previously mentioned, experiments were performed at discrete values 
of reduced frequency to assess the extent to which <a serves as a 
correlating variable for results of the boundary layer velocity surveys. 

As indicated in Table 2, data were obtained at different £ locations for 
reduced frequencies at or near 1.0, 2.2, 3.0, and 4.1. Figures 15 and 16 
show u^/U^ and ^ profiles for boundary layers generated at these values of 
<5. While there are some differences in both the u^/U^ and $ profiles 
within the results at each of the four values of w, it is clear that the 
profiles group more closely than do the profiles at fixed values of 
position alone, Figures 8 and 9, or frequency alone, Figures 13 and 14. 
Thus, within the range of flows investigated, the results tend to group 
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fairly closely when plotted in teems of the y/$ s and at fixed w. 

5.3 Turbulence Level in Oscillating flows 
The measurements made in oscillating flows that yielded the boundary 
layer velocity variation also provided data for assessment of the level of 
longitudinal turbulence. These data form the basis for investigating 
possible time-dependent effects on turbulence in the oscillating flows 
studied. Time-dependent effects can be detected by comparing measured 
levels of turbulence in the boundary layer with those for quasi-steady 
flows, i.e., flows for which no time-dependent effects are present. In a 
quasi-steady (slowly varying) oscillating boundary layer flow, the 
turbulence level at a fixed value of f, and a given y value, as indicated by 
u' r at the point, would vary periodically due to change in the boundary 
layer thickness that results from variation in the free stream velocity. 
This variation can be predicted for a given free stream velocity by use of 
the measured steady flow profiles for u/U m and u' r /u in Figures 4a and 4c. 
From the relation 

5/x - 0.37/(US/v) 0,2 

for an incompressible zero pressure gradient turbulent boundary layer with 

-0 2 

a 1/7 power velocity profile, 5 is nroportional to U . Thus, the 

variation in 8 is 180 degrees out of phase with the free stream velocity 

variation, i.e., an increasing free stream velocity U produces a decreasing 

5. From Figure 4c, it is seen that since y is fixed, this results in a 

variation in u' r /u that is also 180 degrees out of phase with the velocity 

variation. The variation in u, the local boundary layer velocity, is in 

phase with the free stream velocity, as is the variation in u/U. The 

latter can be deduced from Figure 4a and the fac', that the variation in & 

is out of phase with the free stream velocity. In computing u' a (u' /u)u 

* 1 


-25- 


for the varying free stream velocity, u dominates, with the result that the 
variation in u' r is in phase with the free stream velocity. Thus, a 
sinusoidal variation in the free stream velocity in the form U • U m + U^Cos 
b)t produces a quasi-steady boundary layer velocity variation of the form 
u(y,t) » u m (y) + u 1 (y)Cos cat and an essentially sinusoidal variation in u' r 
of the form 

u' r (y,t) - u' rm (y) + u' rl Cos cat (5.1) 

Figure 17 shows the curve of u' vs y/8 computed for steady flow 

l in s 

conditions at 5 * 1.58 m using the profiles for u/U s and u' r /u from Figures 

4a and 4c. Thus the curve represents the measured variation in u' rm with 

y/8 and is applicable to both steady and quasi-steady flows. The 
s 

amplitude u' r ^ depends directly on the amplitude of the free stream 
velocity Uj, which in turn depends on the wedge stroke. Experimental 
results for u' r obtained at three frequencies of oscillation and two wedge 
strokes at £ ■ 1.58 m (results from data sets Dl, D4, D5, D6, and D7) are 
presented in Figure 17 and will be discussed subsequently. The broken 
lines in Figure 17 above and below the u' rm curves show the extremes of 
equation (5.i) determined for free stream conditions corresponding to each 
of the five data sets noted. Data sets Dl, D4, and D6 were obtained with a 
wedge stroke of 1.91 cm, while data sets D5 and D7 were obtained at shorter 
strokes. The quasi-steady amplitude u' r ^ for data set D5 is smaller than 
u' r ^ for data set D4 by approximately the ratio of the amplitudes of the 
free stream velocities, which is given by the ratio of the values for U^/U m 
for the two data sets. From Table 2 the value of this ratio is 0.071/0.121 
= 0.59. Similarly, the quasi-steady value for u' rl for data set D7 is one 
half that for data set D6 because for D7 was one-half that for D6. It 
is useful to compare the amplitude of oscillation of u' to that of the 
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local velocity variation £or quasi-steady oscillating flow. Due to the 
fact that U^/U m for the oscillating flows investigated at £ - 1.58 m was 
relatively small (0.136 or smaller), the ratio u'^/u^ for quasi-steady 
flows at a fixed y value is essentially independent of wedge stroke. 

Figure 18 shows the curve computed for u'^/u^ for quasi-steady flow at £ - 
1.58 m. 

As noted above, an essentially sinusoidal variation was predicted for 

u' r in quasi-steady flows. A sinusoidal variation for u' r is also 

indicated by the measurements of u' r . This is illustrated in Figure 19 in 

which experimental results for u' r and the corresponding fitted cosine 

curves are shown at two values of y/S g for oscillating flows from data set 

D6. Although the results were determined from ensemble-averaging hot wire 

measurements for 210 cycles, there is scatter in the data points about the 

fitted curves. This scatter is directly related to the turbulence level in 

the flow. However, the variations are clearly sinusoidal in nature. It is 

also evident that the variations for u' r shown are not in phase. Fitted 

curves like those in Figure 19 provide a means of displaying the behavior 

of u'“ measured in the boundary layer. It is useful to make comparisons 

with the local boundary layer velocity 

u(y,t) = u m (y) + u 1 (y)Cos[wt + ♦(y)] (1.2) 

by writing for the experimental variation in u' r 

u' r (y,t) = u' rm + u' rl Cos[wt + <Ky) - (3(y) ] (5.2) 

Thus, (3(y) is the phase angle by which u' r lags u. (In quasi-steady 

oscillating flows both 4 and 0 are zero.) In Figure 17, the symbols denote 

measured values of u' „ and the vertical solid lines indicate the extremes 

rm 

of the fitted cosine curves, such as those in Figure 19, as given by 
equation (5.2). For each of the oscillating flows considered, the 
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experlraental results for u' rm differ somevhat from u' rm for quasi-steady 

flow (that measured for steady flow). Except for data set D5> the results 

indicate measured u' values that are slightly above the quasi-steady 

rin 

curve for u' _. Also, for most of the range of y/8 , the amplitudes of 

£ fit 5 

oscillation, u' rl , tend for each data set to differ from those for quasi 

steady oscillating flow. This is particularly evident near y/6 - 1. This 

is also observed in Figure 18 which shows values for the ratio u 'rl /u l 

determined from the measurements for each of the oscillating flows 

considered in Figure 17. The results in Figure 18 are grouped in terms of 

frequency. Although there is scatter in the results at each frequency, the 

data exhibit variations that clearly differ from the quasi-steady 

distribution of u'^/Uj vs y/S g . At 5.04 Hz, the difference is most 

evident near y/8 - 1. The difference between quasi-steady and 
s 

experimental results at the two other frequencies are even more pronounced, 

with the largest values for u'^/u^ occurring at small values of y/8 g and 

at values of y/S near unity. The experimental results for 15.5 and 21 Hz 
s 

show essentially the same trend, and within the scatter of the results, no 
effoct of amplitude of velocity oscillation on the results is evident, as 
indicated by the essential agreement between results for data sets D4 and 
D5 at 15.5 Hz and data sets D6 and D7 at 21 Hz. This indicates that at a 
given point in an oscillating flow u 'rl depends directly on u^. However, 
the ratio u'^/u^ is frequency dependent. 

The variation in the phase angle 8 with y/S gives further evidence of 

S 

time-dependent effects on turbulence in the oscillating flows investigated. 
The measured variations of 8 are shown in Figure 20 for data sets Dl, DA, 
D5, D6, and D7 and are grouped in terms of frequency. The results indicate 
that u' r lags the local boundary layer velocity u by an increasingly larger 
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angle as y/S increases and that large changes in 0 occur across the 

o 

boundary layer. Different variations of |3 with y/8 are indicated for each 
of the three values of frequency. However, there is no evident effect of 
amplitude of velocity oscillation on as indicated by the reasonable 
agreement of results for data sets D4 and D5 as well as for data sets D6 
and D7. It is evident from the collective results in Figures 17 to 20 that 
there are significant effects of unsteadiness on turbulence in the 
oscillating flows investigated at Z * 1.58 m. The time-dependent 
turbulence can be described with reasonable accuracy by use of equation 
(5.2) and the results in the figures. Turbulence levels measured at Z - 
0.61 and 1.12 m exhibit time-dependent effects that differ somewhat from 
those at Z - 1.58 m but are similar in nature. 


6. NUMERICAL METHOD FOR PREDICTING BOUNDARY LAYER BEHAVIOR 


The governing equations for the flows under consideration were taken 
as the ensemble averaged incompressible unsteady boundary layer equations 
incorporating a scalar eddy viscosity relation between stress and strain. 
These equations are 


3u ^ 9v n 

? - + ? - = 0 (6.1) 

3u „ 3u „ 3u 1 3p 3 e) 3u 

?t K 3y""p5x 3y l 1 * vj 5y < 6 ' 2 > 

Following McCroskey and Phillipe (1975) the following coordinate 

transformation is introduced: (x,y, t)-*( 5, t), t) where Z is the distance from 

the turbulent boundary layer origin (the trip wire in Figure 1.), t) = 

(U/25V) 1 , and t = wt. This yields 
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and £ is the dimensionless eddy viscosity. The latter parameter is taken 
in the present study from either the well-known Cebeci-Smith algebraic 
model (1968) or the one-equation model of Glushko (1965). The Glushko 
model uses a single differential equation to describe the transport of 
turbulence kinetic energy in conjunction with an algebraic scale equation 
and includes a time-dependence in the turbulence field as well as the 
velocity field. It is the simplest available model that explicitly admits 
time-dependence in the turbulence field as distinguished from the mean 
velocity field. 

The pressure gradient in equation (6.2) is related to the free stream 
velocity U by 
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and to 8. and 6 in equation (6.4) by 
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( 6 . 6 ) 


Thus it is evident that the solution of the governing equations depends on 
the manner in U varies with position and time in the test section. This 
variation in turn depends on how the flow is generated and may differ 
between experimental facilities. As previously noted, U = U(x,t) can be 
predicted analytically for the present facility as a m M(x,t), where M(x,t) 
is given by equation (2.2). Thus the terms and 8 in equations (6.3) as 
given in equation (6.4) can be determined. 


It is useful to examine the nature of the pressure gradient generated 
in the present flows since the pressure gradient, particularly its phase 
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relative to the free stream velocity, has a significant effect on the 
boundary layer behavior. In this analysis, pressure gradients generated in 
flows with small amplitudes of velocity will be considered. It is 
convenient to write the pressure gradient term, equation (6.5), in 
dimensionless form as 


3P 


o> 3M 

I "Scot 
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+ M 


3M 

mtx 


(6.7) 


where P - p/(kp m ). Substitution of M(x,t) from equation (2.2) into 
equation (6.7) yields 


“ 3x “ Cos f wt + A(x) ~ r (x)l (6.8) 

where 

G(x) - fl-M m S'(x)) 2 + [J- S(x) + M m S(x)A'(x)J 2 ) 1/2 (6.9) 

m 

and 

s<x)[ r + M nA <x>J 

T(x) * tan" 1 _ r sT - ( - ) (6.10) 

m ' ' 

The angle Y by which the pressure gradient leads the free stream Mach 
number (velocity) variation can be determined by comparing the cosine 
'f’s.u'i'.'.t li.viis in equations (2.2) and (6.8). At a fixed x, cos (cot + A) leads 
-Coss(cot A) by n radians. Also, Cos(cot -t- A - T) lags cos(cot + A) by T 
radians. Thus, Cos (cot + A - T) leads -cos (cot + A) by 

Y(x) - n - T(x) (6.11) 

Figure 21 shows G(x) and T(x) in equation (6.8) at three values of 
frequency. At lew frequency the terms S'(x), A(x), and A'(x) in equations 
(6.9) and (6.10) are small, as can be seen from examination of Figure 6. 
Thus the variation of G and T with x is small, T is close to n/2, and Y is 
essentially n/2. Hence, the pressure gradient ~3P/3x leads the free stream 
velocity by n/2. With increasing frequency, G and F depend increasingly on 
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x, which results in Y increasing to values larger than but near n/2. 

The influence of the phase relation between the pressure gradient and 

the free stream velocity on boundary layer behavior can be addressed by 

considering an analysis developed by Patel (1977). Patel considered a free 

stream velocity variation of the form 

U(x, t) - U m + U x e iwtt " < x/Q)1 (6.12) 

in which Q * KU is called the travelling wave convection velocity and U , 
m m 

U^, and K are constants. The pressure gradient produced by U(x,t) in 
equation (6.12) is 
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(6.13) 


For Q • •, equation (6.13) shows that the pressure gradient leads the 
velocity by n/2 radians. Facilities of the type used by Karlsson (1973) 
and Coustiex, et al (1973). (and. as noted, the present facility) produce 
flows in which the pressure gradient leads the free stream velocity by 
approximately n/2 radians. However, Patel's facility generated flows with 
Q * 0.77 U m , which results in a pressure gradient that lags the velocity by 
n/2 radians. As a result, an oscillating-flow boundary layer behavior 
significantly different from that observed in this study or reported in the 
studies of Karlsson and Cousteix was recorded by Patel. For example, Patel 
observed the maximum value for the boundary layer velocity amplitude ratio 
Uj/Uj. at U s 4 to be 1.55, while in the present study this value is seen to 
be about 1.1 in Figure 15d. In addition, Patel recorded large negative 
phase angles in most of the boundary layer, while in the present study 
phase angles are generally both positive negative (see Figure 16). This 
difference is explained as follows. For Q/U m = 1, equation (6.13) 
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indicates that the pressure gradient is zero. For this case only viscous 
forces act in the boundary layer. Patel's analysis shows that this 
produces a boundary layer velocity variation that generally lags the free 
stream velocity and has relatively small amplitude of oscillation. When, 
as in Patel's experimental study, a lagging pressure gradient also acts on 
the flow, the viscous and pressure gradient forces act somewhat in phase to 
produce large amplitudes of oscillation in the boundary layer as well as 
large phase lags relative to the free stream flow. However, a leading 
pressure gradient acts out of phase with the viscous forces and produces 
smaller amplitudes of oscillation and positive values of the phase angle <f> 
except possibly near the boundary layer edge, as observed in the present 
study. 

The numerical method employed to solve equations (6.3) was developed 
by Murphy and Prenter (1981). It uses splined cubic Hermite polynomials to 
represent the stream-normal variation of the flow parameters, with the 
coefficients of these polynomials determined by orthogonal collocation. 

The stream wise derivatives are represented by classical second-order 
accurate finite difference approximations and the temporal variation cum 
linearization is accomplished by either a Crank-Nicolson method or by 
second, order accurate finite difference with Newton-Raphson iteration. The 
choice of iterative or noniterative time marching is made by the user. The 
result is a hybrid finite-element finite-difference scheme which if fourth 
order accurate in y and second-order accurate in x and t. The high 
accuracy in the stream normal direction is consistent with the boundary 
layer assumptions and permits accurate solutions on a relatively coarse 
mesh. 
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7. COMPARISONS OF EXPERIMENTAL AND PREDICTED BOUNDARY LAYER BEHAVIOR 

Computations using both the Cebeci-Smith and Glushko turbulence models 
were performed to predict the behavior for several the oscillating 
turbulent boundary layer cases studied experimentally at £ » 1.58 m. 

Figure 22a shovs the mean velocity profiles predicted by numerical 
computations using both turbulence models. The numerical method with 
either turbulence model predicts slightly higher velocity values in the 
boundary layer than were experimentally observed. Figures 22b and 22c 
compare predicted and experimental results for u^/U^ and the phase angle <p 
at four frequencies. Predicted results are shown for both turbulence 
models. Vhile the predicted and experimental results are not in exact 
agreement for either turbulence model, the results obtained using the 
Glushko model are generally inferior to those obtained from the Cebeci 
Smith model. Therefore the discussion will focus on the results predicted 
using the Cebeci-Smith model. In all cases these results do not show a 
significant influence of unsteadiness beyond y/S * 1, while the 

S 

experimental results do. At frequencies 5.04 Hz and 11.3 Hz the 
experimental and predicted results for both u^/U^ and are in fair 
agreement. At frequencies 15.5 Hz and 21.0 Hz larger differences between 
the predicted and experimental results exist, with the exception of u^/U^ 
at f ■ 21 Hz. Vith the exception of the results at 21 Hz, the predicted 
curves for 4 exhibit increasing values with decreasing y/S for the full 

S 

range of y/S , while the experimental corresponding results (at least the 
results a 5.04 Hz) exhibit a maximum value of <|> at y/S g > 0 and then show 
decreasing values of $ as y/S decreases. The predicted results for <f> at 
21 Hz exhibit characteristics similar to other predicted results for <{>, 
except negative angles are predicted near y/S g = 0.8. For all frequencies 


significant differences exist between the Cebeci-Smith predicted results 
and the experimental results for y/6 less than about 0.1. The behavior of 
the phase angle near the wall is similar to that found by Cousteix (1981). 
Thus there appears to be no reason to question the present results. Near 
the wall the results obtained using the Glushko turbulence model differ 
only in magnitude from those obtained using the Cebeci-Smith model. It is 
likely that the disagreement between predicted and experimental results 
both near the wall and in the remainder of the boundary layer is due at 
least in part to failure of the turbulence models to adequately describe 
viscous effects in the unsteady flows involved. 

8. SUMMARY 

The oscillating flows studied experimentally cover a wide range of 
frequency and reduced frequency. This study was made possible because the 
facility used to generate the time-dependent flows produces an essentially 
pure sinusoidal variation in free stream velocity at all frequencies. The 
study has yielded a quantitative description of the boundary velocity 
variation for y/S greater than about 0.05 at fixed axial positions as a 

5 

function of frequency, fixed frequencies as a function of position, and at 
fixed values of reduced frequency w. Mean velocity profiles, u m /U m vs 
y/5 , were found to exhibit only small deviations from corresponding steady 
flow profiles. The boundary layer velocity amplitude ratio u^/U^ and phase 
angle $ relative to the free stream flow were generally found to depend 
strongly on position and frequency. When plotted in terms of y/S at fixed 
io profiles for u^/Uj^ and <|> tended to collapse toward single curves, 
indicating usefulness of w as a similarity parameter for correlating 
results for the oscillating flows examined. 

The longitudinal turbulence level, as indicated by the measured 


variation in u' r for the ensemble averaged cycle, was found to be strongly 

time-dependent. Quantitative results for turbulence level are presented 

t or flows at l ■ 1.58 m. The mean value of u' r was found to be close to 

that measured for steady flow for all frequencies. However, the amplitude 

of u' and the phase relation between u' and the local boundary layer 
r r 

velocity u exhibited significant differences from those predicted for quasi- 

steady flows. The amplitude of u' was found to vary with y/8 . In 

r 5 

addition, at frequencies above about 15 Hz large amplitudes of oscillation 

in u' occurred both near the wall and near y/8 * 1. The phase relation 

V s 

between u' r and u varied in a pronounced manner across the boundary layer. 

The experimental study of the effects of amplitude of free stream 
velocity on boundary layer behavior was conducted for a limited range of 
V U m . Boundary layer velocity profiles quantities, Uj/U^ and ♦, exhibited 
at most a small influence of amplitude of oscillation. Similarly, u' r /u^ 
as well as the phase angle between u' r and the local velocity in the 
boundary layer showed no influence of amplitude of free stream velocity. 

The numerical method described for predicting boundary layer behavior 
in the time dependent flows studied was applied using the turbulence models 
of Cebeci and Smith, and of Glushko. The required boundary condition at 
the boundary layer edge, the free stream velocity, was supplied using 
analytical expressions developed for the present facility. These 
expressions predict variations in amplitude with position and frequency and 
corresponding phase relations for the free stream flow. Somewhat different 
velocity variations in the boundary layer are predicted by the numerical 
method when the two turbulence models are used. Generally the Cebeci-Smith 
model produces results that are in best agreement with the experiments. 
Numerical results obtained using the Cebeci-Smith model are in many 


respects qualitatively correct. The differences that exist between the 

numerical predictions and measured results are apparently due at least in 

« 

part to the inadequacy of the turbulence models. 
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Table 1. Description of oscillating flow facility. 


Entrance section 

flow-straightening section: 

Diameter; 0.3 m, Length: 0.5 m 
Honeycomb: 0.5 cm hex cell, 3 cm length 

Screen: 0.3 solidity ratio 
Contraction section: 

Length: 0.3 m, Area ratio: 12 to 1 

Test section 

Length: 2.08 m, Cross section: 7.5 cm x 7.5 cm 
Vails: 2.5 cm- thick aluminum 
Boundary layer trips: 0.05 mm dia. wire on four 
walls at x ■ 0.13 m 6 

Unit Reynolds number: 1.52xl0°/m 
Wave generator 

Nozzle throat cross section: 2.13 cm x 7.5 cm 
Wedge: 14 degree included angle 
Wedge strokes: 0.953, 1.27, 1.91 cm 
Frequency: 0 to 35 Hz 


Table 2. Compilation of oscillating flow boundary layer surveys 


Wedge 


Data Set 

5, m 

f, Hz 

b> 

stroke, cm 

V U m 

A1 

0.61 

6.79 

0.54 

1.91 

0.131 

A2 

0.61 

13.1 

1.05 

1.91 

0.142 

A3 

0.61 

15.5 

1.24 

1.91 

0.147 

A4 

0.61 

20.9 

1.67 

1.91 

0.169 

A5 

0.61 

20.9 

1.67 

1.27 

0.112 

B1 

0.91 

19.1 

2.24 

1.91 

0.132 

Cl 

1.12 

3.85 

0.54 

1.91 

0.131 

C2 

1.12 

7.10 

1.00 

1.91 

0.126 

C3 

1.12 

15.5 

2.19 

1.91 

0.138 

C4 

1.12 

21.3 

3.01 

1.91 

0.152 

C5 

1.12 

21.3 

3.01 

1.27 

0.097 

C6 

1.12 

29.4 

4.16 

1.91 

0.239 

D1 

1.58 

5.04 

0,97 

1.91 

0.126 

D2 

1.58 

6.79 

1.31 

1.91 

0.126 

D3 

1.58 

11.3 

2.18 

1.91 

0.118 

D4 

1.58 

15.5 

2.99 

1.91 

0.121 

D5 

1.58 

15.1 

2.91 

1.27 

0.071 

D6 

1.58 

21.0 

4.05 

1.91 

0.136 

D7 

1.58 

21.0 

4.05 

0.953 

0.065 
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Figure 1. Oscillating flow facility. 
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Figure 2. Data acquisition system. 
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Figure 4. Boundary layer velocity and turbulence intensity profiles for 

steady flows: A, 5 = 0.61 m; O, £ = 1.12 m; 0,5= 1.58 m 
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Figure 10. Estimated uncertainties in amplitude of oscillation and phase 
angle for velocity in oscillating flows. 210 cycles. Wedge 
stroke = 1.91 cm. 
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Figure 15. Experimental profiles for u./lL at four values of reduced 
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Figure 19. Experimental variations in u' from data set D6. £ = 1.58 m 

*, ensemble averaged velocity^ 210 cycles; fitted cosine 



£ - 21.0 Hz 


£ = 15.5 Hz 



iriation in (3, the angle by which u' r lags the local boundary 
tyer velocity u, with y/S . K = 1.58 m. 

5 

:periment: O, data set Dl, f = 5.04 Hz, Uj/U = 0.126; A, data 

st D4, f = 15.5 Hz, Uj/U n = 0.121; Q, data set D5, f = 15.5 Hz, 

/U m = 0.071; □, data set D6, f = 21.0 Hz, Uj/U = 0.136; 

i, data set D7, f = 21.0 Hz, U./U = 0.065. 

1 m 




Figure 21. G and F, equation (6.8), vs x at three values of frequency 
for flows in the present facility. Wedge stroke = 1.91 cm. 





